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transport between any two points of such a spacetime. The integrals in this 
formula are evaluated in closed form for a particular spatially-flat model. 
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1. Introduction 

In a recent book [1] , Prugovecki describes in detail a framework for the unification of 
general relativity and geometro-stochastic quantum theory (cf. [2] for a review). The the- 
oretical basis for geometro-stochastic quantization differs from the orthodox quantization 
procedure in that it incorporates the concept of a fundamental length in nature, which 
limits the accuracy with which one can measure the spacetime location of an event. This 
assumption resolves many of the difficulties with divergences encountered in conventional 
quantum field theory. 

In this paper, we focus on the propagation of quantum wave functions for a massive 
spinless boson in a Lorentzian manifold M. In particular, we derive the equations governing 
such propagation in the case that this manifold is a classical Robertson- Walker model of 
our universe. The methods described could easily be applied to other models that have 
received attention in [3] in the context of conventional field theory in curved spacetimes. 

Throughout this paper, we assume that M is space-time orientable, so that it can be 
covered by a global coordinate system of orthonormal frames in which the timelike vector 
is future-pointing and the spacelike vectors form a right-handed triad. In addition, there 
is assumed to be a global time variable t so that we can foliate M into a set of Cauchy 
spacelike hyper surfaces a t - We adopt a metric g with the signature (—1,1,1,1) commonly 
used in relativity theory, and use Planck natural units where c = G = % = 1. 

In section 2 we provide a brief overview of the special-relativistic case in which a 
single wave function on phase space is used to describe the state of the boson. This is 
then extended in section 3 to the case of a globally hyperbolic spacetime through the 
introduction of the general quantum bundle, which is associated to the Poincare frame 
bundle over that spacetime. Parallel transport on this bundle is defined by making use of 
the standard extension of the Levi-Civita connection used in classical general relativity, and 
serves as the basis for the development of the quantum-geometric propagator. In section 
4, we compute the geodesies for the Robertson- Walker model and use them to develop an 
integral representation of the propagators for parallel transport. Finally, in section 5 we 
focus on a particular spatially-flat model for which these propagators turn out to have a 
closed-form expression. 

2. Massive boson states in special relativistic phase space 

Consider a spinless particle of rest mass m > which moves freely in the Minkowski 
space M 4 . The momentum representation of such a particle consists of the space L 2 (V+) 
of the standard wave functions on the forward mass hyperboloid corresponding to m, whose 
inner product is expressed using the usual invariant measure on V+: 



v+ = {keM\ \k° 

As is well-known [4] , the standard interpretation of (p is that of a probability amplitude for 
particle 4-momentum measurements. However, a corresponding relativistically invariant 
probability amplitude does not exist for spacetime measurements [5,6]. 



= / #(A0£2(A0dO«(A0, ^ 2 GL 2 (V+), 
J\r± 



= vV + m 2 }, dn m (k) = 



dk 

2\/k 2 4- m 2 



, k = (k°,k). 



(2.1) 
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Already in the 1930's, various physicists, of whom the most prominent were Heisenberg 
[7] and Born [8], had conjectured that there is a fundamental length I in nature, such 
that no spacetime measurement of an event can be made with greater accuracy than I. 
Measurement-theoretical arguments [9, 10] suggest that this length is in fact Planck's 
length (which is 1 in the natural system of units adopted in this paper). 

In the quantum-geometric framework, this hypothesis forms an integral part of the 
structure of the single-boson Hilbert space F which represents the state of the system. 
In the remainder of this section we present an overview of this theory; for details consult 
sections 3.7 and 3.8 of [1] and chapter 2 of [11]. 

To construct F, we first introduce the concept of the momentum space resolution 
generator, denoted by fj. From a mathematical viewpoint, this resolution generator can be 
chosen to be any wave function from L 2 (V+) satisfying the following two properties; first 
that {fj \ fj) = m/47r 3 , and second that fj is invariant under spatial rotations. 

To understand the physical significance of fj, it is necessary to consider in detail the 
actual process of measuring the spacetime coordinates of a given event relative to some 
inertial frame. Classically, we can imagine that this process is performed with the aid of an 
apparatus employing pointlike particles and light signals. However, a closer examination 
of this procedure shows that the recoil effects produced by the emission and absorption of 
photons will introduce an inherent uncertainty factor. If, however, we replace the pointlike 
particles by extended particles whose momentum space wave function is given by fj, then it 
is possible to define this process in a physically consistent manner which takes advantage 
of the various methods developed in the phase space approach to quantum mechanics [12]. 

The introduction of quantum frames of reference induces a fundamental change in the 
notion of spacetime itself. Specifically, the position and momentum coordinates q G R 3 
and p G R 3 of the particle which are measured at some given time will now be stochastic 
variables (cf. section 3.3 of [1]) whose definitions incorporate Gaussian-type distributions 
X q and Xp on R 3 - The standard deviations of these distributions are of the order of £ and 
i~ x respectively (in Planck units). Hence, it is possible to define a wave function on the 
phase space of the particle which provides a joint probability amplitude for measuring the 
particle as simultaneously having a given stochastic position and stochastic momentum, 
without violating the Heisenberg uncertainty principle. This will be done below. First, 
however, we consider what value the momentum space resolution generator should actually 
have. 

In general, there are infinitely many possible choices for fj. However, in this paper we 
will henceforth consider only the special case 

fj(k) = (m 3 Z,, m p eM-tk°), Z e , m = ^ K 2 (2£m), k = (k°, k) G V+, (2.2) 

where K 2 denotes as usual a modified Bessel function. This case corresponds to the 
unique resolution generator whose phase space representative displays reciprocal invariance 
under transposition of configuration and momentum coordinates, and is the ground state 
eigenvector of Born's [13] quantum metric operator. Hence, its presence reflects the fact 
that stochastic spacetime localization is implemented only to the order of i, whereas the 
corresponding stochastic 4-momentum spread is of the optimal order of magnitude that is 
in accordance with the Heisenberg uncertainty principle [12]. 
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The above choice of r\ induces a mapping W v from L 2 (V+) into the set of complex- 
valued functions defined on the relativistic phase space = {(q,p) £ M 4 x V+}, given 
by 

W v :<p^ip, (p(q,p) = (fj q ,p\ <f>), fj q , p (k) =exp(-iq-k)fj(A~ 1 k), (2.3) 

where A v is the Lorentz boost corresponding to the 4- velocity v = p/m. The image under 
Wrj of L 2 (V+) forms a Hilbert space F with inner product 

(<Pi\<P2)= ¥i{q,P)¥2{q,p)dY, m (q,p), ^(PiGF. (2.4) 

The domain of integration is the 6-dimensional hypersurface E+ = a x V+, where a is 
any Cauchy spacelike hypersurface in M 4 . The hypersurface £+ carries the relativistically 
covariant volume element 

dZ m (q,p) =2p l da l (q)cM m (p), (q,p) e £+. (2.5) 

It can be shown [11,12] that the conditions imposed on the resolution generator fj imply 
that the mapping defined in (2.3) is unitary. Consequently, for normalized <p the values 
assumed by the probability measure 

P„(B)= [ \cp(q,p)\ 2 dZ m (q,p) (2.6) 

JB 

can be interpreted as representing relativistically invariant probabilities for stochastic phase 
space localization. 

This phase space representation has several advantages over the conventional formula- 
tion, which is usually based on wave functions defined on M 4 that satisfy the Klein-Gordon 
equation. In addition to providing a joint probability amplitude for stochastic position and 
momentum as explained above, ip simultaneously satisfies the Klein-Gordon equation in 
its q- variables and the Born-Lande equation [14,15] in its p- variables: 

J* ^J-™ 2 ) V(q, P) = 0, (V ^~ ~ ^ <p(q, P) = 0. (2.7) 

Moreover, on the physical side, (p can be used to define a positive-definite, conserved 
probability current, so that a physically-consistent single-particle theory becomes possible. 
This theory can also be extended to handle particles in external electromagnetic fields. For 
example, it is possible to develop a model in which the particle is coupled to a stochastic 
version of the usual 4-potential Aj, for both the Klein-Gordon equation (section 2.10 of 
[11]), and the Dirac equation [16]. The time development of the wave function in such 
models does not give rise to transitions to negative energy states. Hence, as shown in 
detail in [11] and [12], inconsistencies such as the Klein paradox are avoided. 

The phase space is acted upon by the restricted Poincare group ISOo(3, 1). Each 
element of ISOq(3, 1) can be represented by a pair (6, A), where b is a vector in R 4 and 
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A is a 4-by-4 matrix representing a Lorentz transformation. The action of (6, A) on points 
in 7-*+ is given by 

(6, A)(q,p) = (Aq + 6, Ap), (6, A) e ISO (S, 1), (2.8) 

and has the standard composition law (6, A) (&', A') = (b + Kb 1 ', AA'). In turn, this action 
induces a representation U of I5 , O (3, 1) on F. This representation is defined by 

U(b, A) : <p(q,p) » ^(A-'iq - b) , A" 1 ^, (6, A) G /SO (3, 1), (2.9) 

and physically can be interpreted as providing the relationship between two coordinate 
wave functions corresponding to a particular quantum state, as measured by an arbitrary 
pair of inertial observers in M 4 whose relative motion is determined by b and A. 

The resolution generator defined in (2.2) can be used to construct a special relativistic 
phase space propagator K, which is analogous to the Feynman spacetime propagators 
conventionally used in quantum mechanics. To do this, we first introduce the phase space 
resolution generator rj G F, which is defined in terms of (2.2) and (2.3) by rj = W v fj. From 
the rotational invariance of fj it follows that rj itself satisfies a similar property, i.e. that 
£7(0, Aft)r] = rj for any spatial rotation A^ G £00(3, 1). 

We now define the quantum propagator as follows. Let (q',p') and (q",p") be any 
two points in P+. Then the propagator between these points is given in terms of the inner 
product (2.4) on F by 

K(q",p";q',p') = (r] q ",p" \ Vq', P '), Vg", P " = U(q" \ A v »)ri, = U(q',A v/ )rj. (2.10) 

Using (2.10) it can then be shown that K satisfies all the properties conventionally required 
of propagators. In particular, if (q',p') and (q",p") are arbitrary points in then we 
have 

K(q",p";q',p') = K*(q',p';q",p"), (2.11) 

and 

K(q",p";q',p') = [ K(q",p"; q»',p'») K{q"\p'"- q',p') dS m (q"',p"'), (2-12) 

where, as before, £+ is any allowable hypersurface in P+. From (2.12) it follows that K is 
a reproducing kernel for F, so that the probability amplitudes corresponding to arbitrary 
phase points in are related by 

<p(q",p")= [ K(q",p";q',p')ip(q',p')dZ m (q',p'), <p G F. (2.13) 

Although the quantum propagator enjoys all the usual properties of conventional 
propagators, it has the advantage of being an analytic function as opposed to a distribution. 
Indeed, for the specific choice of momentum resolution generator given in (2.2), it turns 
out that the propagator (2.10) can be evaluated in the closed form 

K(q ,p ;q,p) = r , q = q — q , p = p +p, (2.14) 

m 2 Ze m f£,m{q,P) 
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where Ki again denotes a modified Bessel function, and is given by the Poincare 
invariant expression 



ft,m(q,p) = \ q-q q-p «V'V, q,p e M 4 . (2.15) 

V m m 2 

This result enables us to manipulate propagators algebraically, without having to deal with 
any of the difficulties which arise from attempting to multiply distributions. In particular, 
it will be shown in the next section that the special relativistic propagator can be used to 
construct a quantum-geometric propagator for curved spacetimes. 



3. Quantum propagators in globally hyperbolic spacetimes 

In section 2 we constructed a model of quantum mechanics in special relativistic phase 
space. The system in consideration was represented by a single wave function <p defined 
on the phase space P+. In a general curved spacetime M this is no longer possible in a 
relativistically consistent manner. Instead, it is necessary to construct a representation 
based on various principal and associated bundles over M (cf. sections 2.3-2.6 of [1]). 

Within this context, the set of all possible local states of single-boson quantum systems 
will be represented by a bundle E associated to the Poincare frame bundle PM, with its 
typical fibre the Hilbert space F defined in section 2, and with gauge group ISOq(3, 1) 
acting on F in accordance with (2.9). 

This associated bundle is most easily constructed through the use of a G-product, in 
terms of which it can be written as E = FM XqF. This construction implicitly defines a 
projection map it : E — > M which maps each local quantum state in E into its correspond- 
ing base point in M. For any x G M, the fibre n~ 1 (x) above x is associated to the typical 
fibre F by the soldering map [1,17,18] 

a^:7T-\x)^F, <P X U , u = (a,e l )en- 1 (x), (3.1) 



where IT is the canonical projection map from PM to M which identifies the base point 
in M of any Poincare frame from PM. As usual, the local Poincare frame u has been 
decomposed into a displacement vector a e T X M and a local Lorentz frame (e , e l5 e 2 , e 3 ) e 
(T X M) 4 . 

The physical interpretation of these maps is as follows. The state of a particular 
system corresponds to a section SI/ in E, whose value is given by the quantum-geometric 
propagator described below. Above each point x in M is a local quantum state vector ^ x 
which represents the state of the system at x. The coordinate wave function \I>" : — > C 
of this state as measured by an observer at x is determined by applying to the soldering 
map defined in (3.1) with u = (0, e,), where (e^) is the local Lorentz frame relative to which 
the measurements are performed. 

Since the quantum bundle E is associated to PM, we can use the standard extension 
[19,20] to PM of the Levi-Civita connection on LM to construct a connection on E. The 
connection on PM can be specified by means of the usual 1-forms for infinitesimal Poincare 
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transformations, whose values in terms of the chosen Poincare gauge provided by a section 
s of PM are given by [18] 

^4'( x ) =\ xk [s( e l , [ej,e k ]) + g(ej, [e^e,]) - g(e fc , [e^e,])], 2 ^ 
0* (X) = X* + e* (V x a) , X = X l e t G TM. 

By standard results of differential geometry [19], this connection induces a corresponding 
connection on E, which we refer to as the quantum connection. Below, we examine this 
connection from the viewpoint of coordinate wave functions. 

To do this, we first choose a Poincare gauge s. For any smooth curve 7 joining the 
points x' and x" , the connection on PM gives rise to an operator T y (x", x') : Il~ 1 (x') — > 
Il~ 1 (x") for parallel transport. In turn, this defines the transformation g*(x",x') G 
ISO {3,l) given by 

T y (x",x')s(x') =s(x")g*(x",x'), (3.3) 

where the term on the right side incorporates the (right) group action of ISOo(3, 1) on the 
fibre U~ 1 (x") C PM. In view of the fact (cf. chapter 8 of [21]) that this action satisfies 
the compatibility condition r 7 (x", x')(u'g) = (r 7 (x", x')u')g for any u' G II _1 (x / ) and any 
g G ISOo(3, 1), the transformations defined above satisfy an important composition law. 
Namely, if 7 is smoothly extended to any third point x'" G M, then we have the relation 

g*(x>",x>) = g*(x"',x")g;(x",x'). (3.4) 

In terms of the transformation defined in (3.3), the coordinate wave function of any 
local quantum state in 7t~ 1 (x') C E which is parallel transported from x' to x" along 7 is 
then given by 

= U{ 9 ;(x",x')) rj x '\q,p), (q,p) G V+, (3.5) 

where U is given by (2.9). If we now denote the translation and Lorentz transformation 
parts of g*{x",x') in (3.3) as follows, 

g s(x",x') = (b*(x",x>),As(x",x')), (3.6) 
then, according to (2.9), (3.5) can be written in the form 

<:;,'W) = rf^A^x'^xT'iq-b^x^x')), A^VrV). (3.7) 

The connection constructed on E enables us to generalize the special relativistic phase 
space propagators defined in (2.10) to curved spacetimes. Specifically, for any Poincare 
gauge s, and any two points x' and x" joined by a smooth curve 7, we define the propagator 
for parallel transport from x' to x" along 7 in the gauge s to be a function K^(x", x', £') 
from V+xV+^C given by [1,17,18] 

K*(x",(";x\C) = ( Vc „ I U(g;(x'W))vc)- (3-8) 
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In the above equation, we have introduced the composite phase variable ( = (q,p), so that 
the integration measure dT, m (q,p) on F can be written simply as <iE m (£). This notation 
will be employed throughout the rest of this section. 

In the case where the base manifold M is flat, this propagator is independent of the 
choice of 7. It is then not difficult to show [17] that with an appropriate choice of section 
s, the propagator for parallel transport assumes the same role as the quantum propagator 
considered in section 2. Namely, if we take for s the gauge in which a is identically and 
is parallel to the positive x l -axis, then the propagators defined in (2.10) and (3.8) are 
related by 

K s 1 (x",C';x',C)=K( q " + x",p";q' + x',p'), C = (q>,p'), C = (q",p")- (3-9) 

With curved spacetimes the situation is more complicated. In chapter 4 of [1], Pru- 
govecki proposes the idea of a quantum-geometric propagator, denoted by K. This propa- 
gator is constructed in a similar manner as the conventional Feynman propagator, except 
that the space of paths which is integrated over is formed from arcs of geodesies in M in- 
stead of the usual straight line segments. Physically, K is to be interpreted as representing 
the evolution of a single-particle system against a geometro-dynamic spacetime background 
(cf. section 4.7 of [1] for details), so that, in agreement with (2.13), the quantum-geometric 
wave function 

represents the relative probability amplitude for the detection at x" of the stochastic phase 
space value (" (relative to the gauge s) , for a system prepared at x' in the local state \fv . 

To construct this propagator, we first foliate the given spacetime into a family of 
Cauchy hypersurfaces at indexed by the global time parameter t. Choosing two hyper- 
surfaces and a t >> with t" > t' and some positive integer N, we insert the N — 1 
intermediate hypersurfaces at n , where t n = [(N — n)t' + nt")\ for n = 1, 2, . . . , N — 1. 
For any two points x' and x" lying on the surfaces at' and at" respectively, and for any two 
phase space points (' and Q" (relative to the gauge s), the quantum-geometric propagator 
K s (x", C"; x', C') is defined to be the limit [1,17,18] 

KV',C'V,C) = lim^ / I! ^(xn^-O^'Cn; Xn-l,Cn-l)^0rn,Pn). (3.11) 

Here, 7(x n ,x n _i) is the geodesic joining x n -\ G a tn _ 1 to x n G a tn , the measure on the 
hypersurfaces E = a x V+ is given by 

dT>(x n ,p n ) = 2(p n y d<Ti(x n ) dQ(p n ), (3-12) 

and the intermediate phase points are chosen so that they correspond with the point of 
contact between the tangent space and the base manifold: 

Cn = (-a(x n ),p n ), a l (x n )ei = a(x n ), n = 1, 2, . . . , N - 1. (3.13) 
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The points xo and x n correspond to the initial and final points between which propagation 
takes place, i.e. we have xq = x' , xn = x", (o = C) an d Cn = C" ■ Note also that the 
product in (3.11) does not include an integration over a t ", so that only the intermediate 
hypersurfaces are integrated over. 

In the case of Minkowski space, and with the appropriate choice of gauge, it can be 
shown [1,17,18] using (2.12) and (3.9) that this quantity reduces to the propagator for 
parallel transport given in (3.8). However, very little is known in the case where that 
spacetime is curved. Clearly, a first step to understanding this problem is the evaluation 
of the propagators for parallel transport contained in the product in (3.11). It is this 
problem which we now consider. 

Assuming as always that the resolution generator fj is taken to be the standard one 
given in (2.2), it follows from (2.10) that (3.8) can be written as 

K*(x", C"; C) = {U(q", K»)r] \ U(b + Aq', AA v ,)rj), (3.14) 

where (6, A) = g*(x", x'). In order to evaluate this expression, we use the fact [17] that we 
can write the product of Lorentz transformations in the inner product as AA V > = A\ v rAn, 
where Ar is a spatial rotation. Since rj is rotationally invariant (recall section 2), it follows 
from (2.10) and (2.14) that 

C "; ,', O = K ^YY )] - P = P" + A "' " = ""- A «' - < 3 ' 15 > 

m 2 Z ttm fom(q,P) 

where f^ m is given by (2.15). The calculation of the propagator for parallel transport thus 
reduces to the problem of determining b and A. 

As shown in [18], the 1-forms given in (3.2) can be used to compute an explicit 
formula for g^(x", x'). We shall now present that derivation in a form best suited to the 
computations in the subsequent two sections. 

Suppose that x' and x" in M are connected by a smooth curve 7 with ^{t') = x' and 
7(t") = x" . Choose a positive integer n and divide 7 into n pieces by defining the points 

(n- k)t' + kt" 

Xk=l{tk), t k = , k = 0,1,2, ...,n, (3.16) 

n 

so that in particular xq = x' and x n = x" . By induction on (3.4) we have 



g*(x",x')= Yl g*(x k +i,x k ). (3.17) 

k=n — 1 

Applying the composition law for JS'Oo(3, 1) to the right side and using (3.6) we see that 
for any n = 1, 2, 3, . . ., the total translation and total Lorentz transformation from x' to 
x" are respectively given by 



n-l 



fc=0 



' fc+i 
Yl A*(xi+i,xi) 

l=n-l 



b*(x k+1 ,x k ), A= Yl A*(x k+1 ,x k ). (3.18) 



k=n—l 
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In order to compute these expressions, we use the exponential map from iso(3, 1) 
to ISOq(3, 1). Specifically, let At = (t" — t')/n be the increment in the parameter t 
corresponding to (3.16). Then to the first order in At we have 

b*(x k+u x k ) « AtP s (X fc ), A^xk+uxk) ^exp(AtM s (X fc )), (3.19) 

where = 7 / (£/ c ) is the tangent vector to 7 at t = tf~, and the infinitesimal translation 
and Lorentz transformation are defined in terms of the standard basis elements Pi and 
M lj of iso(3, 1), and the 1-forms given in (3.2), by 



P S (X) = -6UX)Pi, M S (X) = — wf,(X)M^', X G TM. 

2 

Substituting (3.19) into (3.18) gives the approximations 



(3.20) 



n-l 



b « At 



fc=0 



fc+1 



]1 exp(AtM s (X*)) 



.Z=n-1 



P s (X fc ), [Jj exp(AtM s (X fc )). (3.21) 



k=n— 1 



In the limit as n — > 00 the approximations will become exact, and we have as our final 
answer 



b = J T exp [J M S ( 7 '( S )) dsj P s (Y(t)) 
A = Texp ^jT M s ( 7 '(t))^ . 



eft, 



(3.22) 
(3.23) 



4. Quantum parallel transport in Robertson- Walker spacetimes 

As shown in section 3, the propagators for parallel transport are completely deter- 
mined by the quantities b and A given in (3.22) and (3.23). An example of the application 
of these formulae is provided by the classical Robertson- Walker model, whose metric can 
be expressed [22] in the form 



-dt 2 +a 2 (t) 



kr 2 



+ r W 



dQ 2 = d$ 2 + sin 2 6 dcjf , 



(4.1) 



where a(t) is a smooth positive function of t. The parameter k can be either — 1, 0, or 
1, corresponding respectively to spatially hyperbolic, spatially flat, and spatially ellip- 
tic universes. In this section, we compute the geodesies in the Robertson- Walker model 
(employing a method suggested by Prugovecki [23]), and use them to arrive at integral 
representations for b and A. 

To begin, it is convenient to define a new variable R by 



arcsinhr, if k = — 1; 
R = { r. if k = 0; 

arcsin r, if k = 1, 



(4.2) 
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which reduces the metric to the form 



ds 2 = -dt 2 +a 2 (t)[dR 2 + r 2 dn 2 ]. (4.3) 

For the calculation of the geodesies, we concentrate on curves whose spatial projections 
follow paths passing through the origin. Clearly, any geodesic can be obtained from one 
of this form by a suitable shift of the origin. Along such lines 9 and cj) are constant, and 
so the distance along a curve R(t) ,t'<t<t" can be expressed as 

1 = £ £(R,R,t)dt, C(R,R,t) = ^±(l-a 2 (t)R 2 ), (4.4) 

where the plus sign is taken for timelike geodesies and the minus sign for spacelike geodesies. 

To minimize or maximize / as a functional of the path we use the Euler-Lagrange 
equation, which in the case of timelike geodesies takes the form [24] 

dC _ d dC _ d I a 2 (t)R 
~dR~ Jt~dR ~ dt y^i_ a 2( t ) R 2* 

Integrating with respect to t gives 

R(t) =C f ; dt (4.6) 

J a{t)^C 2 + a 2 {t) 

where C > is a constant of integration. A similar computation for the spacelike case 
yields 

R(t) =C f ; dt (4.7) 

V ' J a(t)^/C 2 -a 2 (t) V ^ 

We now turn to the problem of calculating b and A. For simplicity, we adopt a 
Poincare gauge in which the displacement vector a is identically 0. From (4.3) it follows 
that a natural choice for the basis vectors is given by 

11 1 

e = d t , ei = -7-rdn, e 2 = — j-rde, e 3 = — — . d^. (4.8) 
a(t) rait) ra[t) smd 

Their Lie brackets are then computed to be 



(4.9) 



r i a '(t) r ! a'(t) r i 

e , ei = 7TT- e l' e 0' e 2 = 77V" e 2' e o, e 3 = 77V e 3' 

a(t) a(t) a(t) 

r . 1 dr r . 1 dr r , cot 9 

ei,e 2 = 7TT-7^ e 2 5 ei,e 3 = 7^T5 e3 > e 2,e 3 = --e 3 . 

ra{t)dR ra{t)dR ra(t) 

From (3.2), the nonzero coefficients corresponding to Lorentz boosts are found to be 

^ 1 (e 1 )=^ 2 (e 2 )=^ 3 (e 3 ) = ^, (4.10) 
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while for spatial rotations they are 



tv s 12 (e 2 ) =(4,(e 3 ) = 



1 dr 
ra(t) dR' 



^2 3 ( e 3) = 



COt# 

r a(t) 



(4.11) 



Since the displacement vector a is assumed to be identically 0, it follows from (3.2) that 
the 1-forms for translation will simply be given by #g(X) = X 1 . 

For the computations of b and A we will take 7 to be one of the special geodesies 
considered in (4.4). Because of the similarity of (4.6) and (4.7), it is convenient in the 
following calculations to consider the cases of timelike and spacelike geodesies together. 
Following the example of (4.4), we will adopt the convention that whenever a ± sign occurs 
the plus refers to the timelike case and the minus to the spacelike case. 

Let the initial and final points of 7 be x' and x", corresponding to the global time 
values t' and t" respectively. From (4.6) and (4.7) the tangent vector is seen to be 



l\t) = e + Jyei, c(t) = y/C*±a*(t). 



(4.12) 



It then follows from (3.20) and (4. 10)- (4. 12) that the infinitesimal translation and Lorentz 
transformation elements of iso(3, 1) along 7 are respectively 



(4.13) 



Clearly the quantities M s (^f'(t)) commute for all values of t, and so (3.22) and (3.23) 
become 



b=- 



exp — 



A = exp 



* Ca'(s) 



M 01 ds 



a(s)c(s) 

Ca'{t) 



C 

1,— ,0,0 



a(t)c(t) 



c(t) 
M tyl dt\. 



dt, 



(4.14) 



(4.15) 



The case a(t) = 1 corresponds to Minkowski space. Here, (4.14) and (4.15) reduce to 
the equations 

^0 + ~t£^Pi 



b=(t'-t") 



VC 2 ±1 



1. 



(4.16) 



In this case it follows from (3.7) that the coordinate wave functions at the endpoints are 
related by 

C 



Po + 



Vc 2 ±1 



Pi 



p 



(4.17) 



confirming the fact that two observers in flat space at rest relative to one another and with 
coincident axes will measure the same wave function for the system, modulo a spacetime 
translation. 
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5. Explicit propagators for a spatially-flat Robertson- Walker model 

The integral representations for b and A which were derived in section 4 generally 
cannot be evaluated in closed form. Indeed, even in the spatially-flat case corresponding 
to R = r, the determining equations for C which are obtained from (4.6) or (4.7) cannot 
be solved explicitly, except in special cases. In this section we consider one such case, the 
rapidly-expanding universe corresponding to a(t) = t. 

We begin by choosing coordinates (t, 3 ) for M, in terms of which the metric 

assumes the form ds 2 = —dt 2 + 1 2 ((dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ), and we take as our Poincare 
gauge a section s of FM in which the displacement vector a is identically and the 
orthonormal frame is given by 

d Id Id Id . ^ 

The initial and final points of the geodesic are taken to be x' = (t', x') and x" = (t", x"), 
with t" > t' . We decompose the spatial separation of the two points as x" — x' = pn, where 
p > and n = (ni, 77-2,^3) is a unit vector. For future reference we define the quantities 
Pnull and psi by 

Pnuii = In y , pa = In I ~, I • (5.2) 

Consider first the case where the geodesic joining x' to x" is timelike. We begin by 
calculating the constant of integration in (4.6), which for definitiveness will be denoted as 
C t \. For this purpose it is convenient to define the new variable u = >Jt 2 + C 2 V with u' 
and u" defined accordingly. Using (4.6) with a(t) = t we get 

C - S = ^(V + C «\f- C J<l). (5.3) 



It- ts/Cl + P 2 V(«'-Cu)(u" + C,i). 

Exponentiating and squaring this equation gives the pair of relations 

(u' + C tl )(u" - C t i) = e p t't", (u' - C t i)(u" + C t i) = e-W, (5.4) 

which when added together lead to the result 

e~p/ 2 (e 2 P -l)t't" , . 

Cti = — V ; < p < p QU n. 5.5 

2^/{t" -en>)(ePt" -t>) ^ ' 

We next proceed to the calculation of A. To simplify this computation, it is convenient 
to first consider the case in which x 2 and x 3 are constant along 7, so that n = (1, 0, 0). In 
this case the Lorentz transformation is found from (4.15) to be 

A = e W [-£ t J^__ M 01 dt\, n= (1,0,0). (5.6) 
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A comparison of this formula with (5.3) shows that A is the standard Lorentz boost in 
the x 1 direction, with coshp in its (0,0) and (1,1) entries and — sinhp in its (0,1) and 
(1,0) entries. This result can then be extended to arbitrary values of n by noting that 
since the spacetime in question is spatially flat, it is only necessary to conjugate A with an 
appropriate spatial rotation matrix which maps (1, 0, 0) into n. In this manner the Lorentz 
transformation for a general geodesic 7 is found to be 



A 



coshp —ni sinhp —772 sinhp —773 sinhp 

—ni sinhp nf(coshp — 1) + 1 ni?72(coshp — 1) ni?i3(coshp — 1) 

—772 sinhp n2ni(cosh p — 1) n|(cosh p — 1) + 1 772773 (cosh p — 1) 

—77,3 sinhp 77,377,1 (cosh p — 1) 77.377.2 (cosh p — 1) n|(cosh p — 1) + 1 . 



(5.7) 



To calculate the translation b provided by (4.14), note that due to the spatial flatness 
of M the spatial part of b will be parallel to the direction of propagation. Thus, if we 
denote the temporal and spatial parts of b by br and 65, then b can be written in the form 
b = (6t, M)- In order to compute these quantities, we again consider first the case of 
propagation along the x 1 -axis, for which n = (1, 0, 0), and obtain from (4.14) the formula 



b = 



exp 




0,0 



dt. 



(5.8) 



Using (5.3) to evaluate the matrix integral inside the exponential, we get by equating the 
first two components of the left and right sides of (5.8) the equations 



bx = 



cosh(lnf) 



C t i sinh(lnv) 



u 



dt, 



C t i cosh(lnf) 



u 



— sinh(lni> ) 



dt, 



Cti 



(5.9) 



These integrals are most easily computed by expressing them in a form in which the limits 
are u' and u" . In the case of 6y, this change of variables produces an elementary integral, 



br = — 



C tl (v 2 -l)-u(v 2 + l) 



U ' 



u , u"(u'-u") 
du = 



(5.10) 



while for 65 a similar calculation gives 



u(v 2 



1) - C tl (v 2 + 1) C t i 

: dU = 



du = 



C t i(u' 



u 



(5.11) 



Substituting u = ^/t 2 + C t 2 L into (5.10) and (5.11) and using (5.5) to eliminate C t \ gives 
a pair of expressions involving e p . These can be simplified by converting the exponentials 
into hyperbolic functions, and the final result is then found to be 



br = t' coshp — t", 



bs = —t' sinhp. 



(5.12) 
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Using the spatial flatness of M, it is not difficult to see that these formulae in fact hold 
for all values of the direction vector n. 

The computations performed above can be repeated for the spacelike case, with similar 
results. However, one important difference arises. To see this, consider first the determi- 
nation of the constant of integration in (4.7), denoted this time by C s \. Substitution of 
a(t) = t into (4.7) gives 

t" 

P = J . 7=^ ^ ( 5 - 13 ) 



A comparison of this equation with (5.3) immediately gives the solution C s \ = iC t i, with 
C t i given by (5.5). Note that C s \ is actually real since in the case of spacelike geodesies 
we have p > p nu \i and hence t" — e p t' < 0. However, it follows also from (5.13) that the 
range of C s \ will be restricted to the interval [i",oo). The value C s \ = t" corresponds to 
the greatest possible spatial separation of the points x' and x", which is achieved by a 
geodesic with points given by {(i, x' + r max (t)n) 1 1' < t < t"}, where 



f- 

Jt' s 



t" 



-ds = — In 



(t» - ^/(t»y - {t'Y) (t» + v/(t") 2 -^ 2 ) 



(5.14) 



Substitution of t = t" into this equation and comparison with the second equation in (5.2) 
reveals that such a geodesic will achieve a spatial separation of magnitude p s \ between x' 
and x" . Consequently, geodesies connecting two points with a spatial separation greater 
that p s \ are formed by smoothly extending the geodesic given by (5.14) into a straight line 
with spatial projection parallel to n, lying in the Cauchy hypersurface £ t » = {(£", x) | x e 
R 3 }. This extension is easily seen to satisfy the geodesic equation at every point. 

Let us first consider the geodesies for which p nu n < p < p s \- In this case the translation 
and Lorentz transformation for geodesies whose spatial projection is parallel to the x 1 axis 
are found from (4.14) and (4.15) to be 



b=- 



exp - 




C s i 



,0,0 



dt, 



(5.15) 



exp - 



M 01 dt 



(5.16) 



Substitution of iC t \ for C s \ in the above expressions and comparison with (5.6) and (5.8) 
shows that the formulae for b and A in the spacelike case formally reduce to those for the 
timelike case. 

For geodesies 7 corresponding to p > p s \ we define an intermediate point x = (i, x) on 
7 by i = t" and x = x' + p s \ii, so that as x varies from x' to x, 7 follows the geodesic given 
by (5.14), while as x varies from x to x", 7 follows a straight line along the hypersurface 
Yit>>- It then follows from the composition law for 7S'Oo(3, 1) that the translation and 
Lorentz transformation from x' to x" are given by 



A = A*Jx", x)A s Jx, x'), b = b s Jx", x) + A s Jx", x)b s Jx, x'). 



(5.17) 
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The quantity A*(x",x) can be obtained from (3.23) by considering t to represent spatial 
instead of temporal distance, and integrating from p s i to p along the tangent vector t"e\. 
An easy computation then shows that A*(x",x) is a matrix of the form (5.7) with coshp 
and sinhp replaced by cosh(p — p s \) and sinh(p — p s \) respectively. Consequently, it follows 
from the first equation in (5.17) that the total Lorentz transformation from x' to x" is 
given as before by (5.7). 

To calculate the total translation, we decompose b^(x, x') and b^(x", x) into temporal 
and spatial parts in the earlier-mentioned manner. In terms of this decomposition, the 
second equation in (5.17) can be expressed in the matrix form 



bs 



brix", x) 
b s (x",x) 



+ 



cosh(p - psi) - sinh(p - p s \) 
- sinh(p - p s i) cosh(/9 - p s i) 



br(x, x') 
b s (x,x') 



(5.18) 



To calculate brix.x') and bs(x, x'), first note that from (5.2) we have 



cosh/9 s i = t —, 



sinh psi 



Substituting p = p s \ into (5.12) and applying the first formula above then gives 

brix^x') = 0, bs(x, x') = — t'sinhp s i. 



(5.19) 



(5.20) 



Using (3.22), the quantities bT(x",x) and bs(x",x) can be evaluated by a calculation 
similar to the derivation of A^(x", x) as 



b T (x", x) = t"[cosh(p - p 8 i) - 1], b s (x", x) = -t" sinh(p - p sl ), 



(5.21) 



and so we get by substituting (5.20) and (5.21) into (5.18) that for p > p nu u the spacetime 
translation is given by 



b T = t"[cosh(p - p s i) - 1] + t' sinh(p - p a \) sinhp s i, 
bs = -t" sinh(p - psi) - t' cosh(p - p s i) sinhp s i. 



(5.22) 



Expanding out the hyperbolic functions in these expressions and then using (5.19), we find 
that once again they reduce formally to (5.12). As before, this result is easily seen to hold 
for all choices of the direction vector n. 

So far, we have not made any mention of the case where the geodesic joining x' to x" 
is null. Although this case could be handled directly using (3.22) and (3.23), it is more 
straightforward to note that since b and A are continuous functionals of the path 7, these 
quantities will be given by the formulae already derived. 

In conclusion, the expressions in (5.7) and (5.12) provide, through (2.15) and (3.15), 
an explicit form of the propagator for parallel transport between any two points x' and 
x" in M for which t" > t' . In the concluding paragraph of [18], Prugovecki speculates 
that the quantum-geometric propagators associated with rapidly expanding or contracting 
universes may lead to violations of strict Einstein causality. The calculations presented 
above can serve as a basis for the investigation of such phenomena. 
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